Boundary susceptibility in the spin- 1/2 chain: 
Curie like behavior without magnetic impurities 



o 



(N 



C3 

s 

I 



O 
O 



O 



X 



Satoshi Fujimoto 1 and Sebastian Eggert 2 
1 Department of Physics, Kyoto University, Kyoto 606-8502, Japan 
2 Institute of Theoretical Physics, Chalmers University of Technology, S-41296 Goteborg, Sweden 

(February 2, 2008) 



We investigate the low-temperature thermodynamics of the spin- 1/2 Heisenberg chain with open 
ends. On the basis of boundary conformal field theory arguments and numerical density matrix 
, renormalization group calculations, it is established that in the isotropic case the impurity suscep- 

tibility exhibits a Curie-like divergent behavior as the temperature decreases, even in the absence 
of magnetic impurities. A similar singular temperature dependence is also found in the boundary 
contributions of the specific heat coefficient. In the anisotropic case, for 1/2 < A < 1, these bound- 
ary quantities still show singular temperature dependence obeying a power law with an anomalous 
dimension. Experimental consequences will be discussed. 



Impurity doping in low dimensional antiferromagnctic systems has been a topic of great interest in recent years 
[1-11]. New interesting boundary phenomena give a better physical insight to the underlying structure of the strongly 
correlated states in those systems. In one-dimensional systems impurities typically cut the chains and play the role 
of effective boundary conditions. The open ends give rise to intriguing boundary phenomena, which result in a 
characteristic temperature dependence of the excess susceptibility and specific heat due to the impurities. Famous 
examples include the spin-1 chain where effective spin-1/2 degrees of freedom are created near the open ends which 
show clear experimental signatures [1]. Even in the two dimensional Heisenberg model it has been postulated that 
non-magnetic impurities may give rise to a divergent Curie- like impurity susceptibility [2]. However, the temperature 
dependence of the boundary susceptibility has been much less clear in the case of the antiferromagnetic spin-1/2 chain, 
which is possibly the oldest and most studied prototype of a strongly correlated system. According to the Bethe ansatz 
solutions for integrable systems with open boundaries, the boundary part of the uniform spin susceptibility at zero 
temperature behaves like ~ l/[/i{ln(/i)} 2 ] for a small magnetic field h in the case of SU(2) spin rotational symmetry 
[3-6]. This divergent behavior implies that in the vicinity of the open edges spin excitations are very sensitive to 
' small external fields. It seems therefore likely that the susceptibility should also show divergent behavior for small 
temperatures at zero field. On the other hand, irrelevant boundary operators are known to only produce a small 
finite impurity susceptibility at low temperatures [7]. We now show that an addition to the surface energy from the 
marginally irrelevant bulk operator actually gives rise to the leading singular temperature behavior of the impurity 
susceptibility and the specific heat coefficient corresponding to a Curie-like behavior with logarithmic corrections. 
We consider the antiferromagnetic spin-1/2 Heisenberg chain with open boundaries 
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H XXZ = JJ2 i S i S f+i + StSf +1 + AS?S? +1 ] (1) 

S3 , 

in the limit N J/T, i.e. the thermodynamic limit. Note that in the opposite limit T/J <C 1/N the behavior is 
O ■ trivially described by the ground state, which is a singlet for even N with exponentially small susceptibility as T — > 
and a doublet for odd N with a Curie law behavior. The crossover to the ground state behavior is experimentally 
important and will be discussed later, but for now we will consider the low temperature behavior in the thermodynamic 
limit 1/N < T/J < 1. In the massless region, < A < 1, the low energy fixed point of the Hamiltonian (1) is the 
■ Tomonaga-Luttinger liquid, which belongs to the universality class of the Gaussian theory with the central charge 
c = 1. The low energy effective Hamiltonian with leading irrelevant interactions has been exactly obtained by 
Lukyanov [12]. Using this effective theory we can implement a perturbative expansion of the free energy in terms 
of leading irrelevant interactions, and evaluate impurity corrections of order 1/N. Following the idea of Cardy and 
Lewellen [13], we consider a semi- infinite cylinder on which the direction tangent to the circumference is taken as an 
imaginary time axis, and the direction perpendicular to it as a space axis. The circumference is equal to the inverse 
temperature 1/T. We define the boson phase field on this geometry, of which the mode expansion is given by 



c (x,r) = Q-inTPx + 2wT^L= + 1 V ha n e-^ T < x+i ^ + a n e~ 2 * Tn ^- iT \ (2) 
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where w is the winding number of <jf , K is the Luttingcr liquid parameter, and the operators satisfy the commutation 
relations [a n , a m ] — [a n , a m ] — n5 n + m .o, [a n , a m ] — 0, and [Q, P] = i. Then, the low-energy effective Hamiltonian on 
the semi-infinite cylinder is written as, 

H C = H^ + H^ nt , (3) 

H c = f o 1T ^[(dr<t> c ) 2 + (d x <P c ) 2 ), (4) 
^ T dr 

■ 
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H fnt = ~a 2K ~ 2 \ I —cos(V8K<P c ). (5) 







Here the constants K and a are parametrized as, K = [1 — (1/ir) cos _1 (A)] _1 , a = 2(K — I) / [J K s\n{-n / K )], and A is 
given by equation (2.24) in Ref. [12]. For this choice of the lattice parameter a, the velocity of spinons is unity in the 
unit of 1/a. Thus the Hamiltonian (3) is obtained by simply interchanging time and space coordinates of the usual 
Hamiltonian defined on a circumference of L — aN. 

We express the partition function by using the transfer matrix exp(—LH c ) and the boundary state for H$, \B). 
The lowest order terms of the free energy for the Hamiltonian (3) are given by, 

F = _^l n( 0|e- i ^| J B ) + ^ f L dx (Q\eM-LH§)eMxHS)Hf nt exp(-xHS)\B) | ^ (g) 
L L J {0\exp(—LH§)\B) 

where |0) is the ground state of H§. The first term on the right-hand side of Eq. (6) is the free energy of the c = 1 
Gaussian model. The second term is the 1/L correction that emerges as a result of boundary effects. For the periodic 
boundary condition, this term vanishes. In the case of an external magnetic field h, the free energy is evaluated by 
shifting the boson field <j) c (x) to (j) c {x) — <j) c (x) — y/K/2hx. Using Cardy and Lewellen's method [13], we compute the 
first order term as, 

A (iraT) 2K (&\B) f L cos(2Khx) 



L 2ira (0\B) J [sinh(27rTa;)] 

where |$) is the primary state that corresponds to the conformal field exp(i\/8K(j) c ). In the following, we evaluate 
Eq. (7) exactly for the free open boundary condition. 

To calculate the prcfactor in Eq. (7), we utilize properties of the boundary state. A conformally invariant boundary 
condition is imposed by demanding T(z) = T(z) at the boundary [14]. Here T(z) (T(z)) is the holomorophic (anti- 
holomorophic) part of the stress energy tensor. For the Gaussian model with c = 1, this condition leads to the 
following constraint on the boundary state [15-17], 

(a n ±a- n )\B) =0, (8) 

where the plus (minus) sign corresponds to the Neumann (Dirichlet) boundary condition. Equation (8) is solved in 
terms of the Ishibashi states, which, in our case, are constructed from the highest weight states of the U(l) Kac-Moody 
algebra \v,w) and their descendants [18]. Here v, w are integers specifying the U(l) highest weight state. The leading 
irrelevant interaction (5) is expressed by the primary field exp(iV&K(j) c ) which corresponds to the primary state |2, 0). 
Therefore, is non-vanishing, only if \B) contains |2,0). The Neumann boundary state which consists of the 

highest weight state |0,tu) and its descendants does not satisfy this condition. On the other hand, the Dirichlet 
boundary state, 

/ K\ 1 / 4 00 — 

^ e-^e-S-i 1 "-^-^^^), (9) 

has a finite overlap with |<f>). Then, we have, 

v ' 1 ' = l (10) 
(0\D) [ Uj 

Carrying out the integral in Eq. (7) and using Eq. (10), we obtain corrections to the boundary part of the free energy, 

^ B = _ 2^N {27TaT)2K ~ lRc[B(K + 'Of' 1 _ 2K)l 
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where B(x, y) — T(x)T(y)/T(x + y). 

The boundary contribution to the spin susceptibility is derived from Eq. (11), 

XB = ^Sf B{K ' 1 - 2K)[7r2 - w( K m™Tf K -\ (12) 

with tp'(x) = dip(x)/dx. Note that for 1 < K < 3/2 (1/2 < A < 1), the boundary spin susceptibility xb shows a 
divergent behavior <~ 1/T 3 ~ 2K , as temperature decreases. This anomalous temperature dependence is also observed 
in the boundary part of the specific heat coefficient computed as, 

°E = ^ {2K 1)(2K 2)B{K ! _ 2K)(2iraTf K - 3 . (13) 

Boundary terms that are regular in h and T give higher order corrections and have been neglected here. We would 
like to stress that in the formulas (12) and (13) there is no free parameter, and the prefactors are exact. The divergent 
behaviors of Eqs. (12) and (13) for T — ► are physically understood as follows. In contrast to the bulk Heisenberg 
chains, the ground state degeneracy at the boundaries gives rise to large spin fluctuations, which disturb the spin 
singlet formation. It should be emphasized that the singular behaviors arc not due to the presence of boundary 
operators, but interpreted as a consequence of finite-temperature corrections of the surface energy and the boundary 
entropy ln(0|B) caused by bulk irrelevant interactions. 

At zero temperature with a small magnetic field, a similar singular behavior appears in the field dependence of the 
boundary spin susceptibility given by, 

Xb(T = 0) = ^ fj^ 1 sin(^)r(3 - 2K)h 2K -\ (14) 

The zero temperature susceptibility is also derived from the Bethe ansatz exact solution by using the Wiener-Hopf 
method. We have checked that Eq. (14) coincides with the result obtained by the Bethe ansatz method. 

Now let us consider the isotropic case K = 1. The free energy correction (11) possesses poles for K = 1. To deal 
with these singularities, we follow the procedure considered by Lukyanov for bulk spin systems [12]. We rewrite Hi nt 
in terms of the SU(2) current operators, 

Hint = J ^[ff|| JoJo + Y (J + J - + J - J +)l- ^ 

The exact expressions for the running coupling constants are known as g^ = 2(1 — 1/K)(1 + q)/(l — q), g± = 
4(1 — l/i^)g 1 / 2 /(l — g) [12,19]. Here the parameter q is the function of T and h, of which the expression is also known 
exactly. On the other hand, for the value of K close to 1, Eq. (11) can be expanded in power series of 1 — 1/K. 
Comparing the expansion of Eq. (11) with the expression for g\\ and g±, we can write the free energy correction (11) 
as a power series expansion in terms of g\\ and g±. Then, taking the limit K — > 1 and g\\ ,g±^>g 7 we end up with, 

for h «T. The running coupling constant is determined from the equation, 

g- 1 + \ Ma) = -Re[v(i + ^)] + Msjl^l (17) 

where ip(x) is the di-gamma function. Using Eqs. (16) and (17), we obtain the leading term of the boundary spin 
susceptibility and the specific heat coefficient, 

g g 2 3V>"(1) 1 / lnln(q/r) \ 

XB ~ VINT + 12~/VT ( 2tt 2 ' + "' ~ l2NT\n(a/T) \ ~ 21n(a/T) + "7' ( ' 

Cb _ 9 2 , 5 5 3 , 1 f 1 lnln(a/T) | \ (1Q) 



T 2NT ANT "' 2NT(ln(a/T)) 2 \ ln(a/T) 

where a = ^/-7r/2cxp(l/4 + 7) with 7 the Euler constant. At zero temperature, the boundary spin susceptibility for 
a small magnetic field is also derived from the opposite limit /i>Tof Eq. (11), 
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Xb(T = 0) = 



5g 3 

ANh ^ 8Nh 
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lnln(27ra//i) 



(20) 



4Nh(\n(2ira/h)) 2 \^ ln(2wa/h) 

This result (20) coincides completely with that obtained by the Bethe ansatz exact solutions [3-6]. 

The prefactor of the leading terms in Eqs. (18), (19), and (20) are exact for the isotropic Heisenberg chain, but 
not universal. For example a frustrating nearest neighbor coupling J2 will change the prefactors and at the critical 
point J2 ~ 0.241167 [20] the singular behavior is absent. Special logarithmic singularities at boundaries have also 
been found in the NMR relaxation rate [8,9], which are also directly related to the bulk marginal operator. However, 
in Eqs. (18) and (19) it is even the leading T-dependence that is changed due to this irrelevant operator. 
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FIG. 1. The boundary susceptibility \b as a function of T according to the DMRG calculations (black points), compared 
to the result of the field theory in Eq. (18). Inset: approximate temperature dependence of the average Curie constant per 
impurity for an impurity density of p = 0.1% 

For comparison we have also used the numerical density matrix renormalization group for transfer matrices (TMRG) 
applied to impurity problems [10]. This method can calculate local expectation values and the impurity free energy 
directly in the thermodynamic limit N — > 00. When evaluating the impurity susceptibility we need to take the second 
derivative of the impurity free energy and therefore subtract two large numbers, which becomes inaccurate for very 
low temperatures. For the lowest temperatures (T < 0.1 J) we have therefore instead summed over the excess local 
responses in a range around the open ends, which gave more accurate results. This method agrees with taking the 
second derivative for higher T and should also be a good approximation for T < 0.1J. Note, however, that taking 
the excess response only at the site closest to the open end as was done in Ref. [10] gives a weaker temperature 
dependence. The results are shown in Fig. 1 without any adjustable parameters. 

Finally, we would like to remark on the implications of our findings for experimental observations. In experimental 
quasi-one dimensional systems such as Sr 2 Cu0 3 a small density p of impurities is always present in the form of intrinsic 
defects of the crystal which effectively cut the chains. A corresponding Curie contribution has been observed that can 
be strongly reduced by careful annealing [21,22], which implies that this Curie tail cannot be due to magnetic impurities 
in the sample. For extremely low temperatures T/J <C 1/N such a Curie behavior can be explained by finite chains 
with odd ./V that have locked into their doublet ground state [11]. We have now shown that a Curie- like behavior can 
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even be expected for higher temperatures, albeit with a different Curie constant in Eq. (18). The crossover between 
the ground state and the thermodynamic behavior is well understood, and in fact the partition function can be 
written explicitly if irrelevant operators are ignored [7] . Roughly, the crossover occurs when the temperature becomes 
comparable to the finite size energy gap T ~ irv/N. For a carefully annealed sample of Sr2Cu03 with p ~ 0.013% 
[21] this means that the ground state contribution is only significant for T < 0.001J <~ 2K. For higher temperatures 
the impurity susceptibility is dominated by the expression in Eq. (18). Experimentally this means that the effective 
Curie constant first drops logarithmically as the temperature is lowered. In previous studies this slight change of 
the Curie constant has been fitted to a Curie Weiss correction [22] assuming a phenomenological interaction between 
the impurity moments. The result in Eq. (18) now gives a clear prescription for the expected form of the impurity 
susceptibility. At the lowest temperatures T < hpJ the Curie constant increases again sharply to a limiting value 
of p/8J as T — ► due to the ground state contributions of the chains with odd N [11], leading to a characteristic 
minimum in the average Curie constant Txavg- The approximate behavior of the averaged impurity susceptibility 
Txavg is shown in the inset of Fig. 1 for p = 0.1% by averaging over all chain lengths and assuming a sharp crossover 
from ground state to thermodynamic behavior. 

In summary, we have studied the boundary thermodynamics for spin- 1/2 Heisenberg chains with open ends by 
using boundary conformal field theory and numerical TMRG calculations. It has been shown that the boundary 
contributions of the spin susceptibility and the specific heat coefficient exhibit divergent behaviors as the temperature 
is lowered, which explains experimental observations in quasi one dimensional compounds. 
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